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Abstract 

^ ■ We calculate p,7r^,K^ and A(+S'^) rapidity distributions and compare 



to experimental data from SIS to SPS energies within the UrQMD and 



^ ■ HSD transport approaches that are both based on string, quark, diquark 

{q, q, qq, qq) and hadronic degrees of freedom. The two transport models do 
not include any explicit phase transition to a quark-gluon plasma (QGP). 
It is found that both approaches agree rather well with each other and with 
the experimental rapidity distributions for protons, A's, vr^ and K^. In- 
spite of this apparent agreement both transport models fail to reproduce the 
maximum in the excitation function for the ratio /tt^ found experimen- 
tally between 11 and 40 A-GeV. A comparison to the various experimental 
data shows that this 'failure' is dominantly due to an insufficient descrip- 
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tion of pion rapidity distributions rather than missing 'strangeness'. The 
modest differences in the transport model results - on the other hand - can 
be attributed to different implementations of string formation and frag- 
mentation, that are not sufficiently controlled by experimental data for the 
'elementary' reactions in vacuum. 
PACS: 24.10.-i; 25.75.-q; 11.30.Rd; 13.60.-r 

Keywords: Nuclear-reaction models and methods; Relativistic heavy-ion collisions; Chiral 
symmetries; Meson production 
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I. INTRODUCTION 



The ultimate goal of relativistic nucleus-nucleus collisions is to reanalyse the early 'big- 
bang' under laboratory conditions and to find the 'smoking gun' for a phase transition from 
the expected initial quark-gluon plasma (QGP) to a phase characterized by an interacting 
hadron gas [1-3]. Though evidence for a 'new phase of hadronic matter' at the SPS has 
been claimed [4] , a direct proof - according to the understanding of the authors - is still 
lacking [5,6]. Furthermore, nucleus-nucleus collisions with initial energies per nucleon of 
Ri 21.3 A-TeV (v^ = 200 GeV) are available now at the Relativistic-Heavy-Ion-CoUider 
(RHIC) in Brookhaven, an order of magnitude higher than at SPS energies {^/s « 17-19 
GeV). In central collisions of Au -|- Au nuclei here energy densities above 4 GeV/fm^ 
are expected [6] . These estimates are based on the Bjorken prescription [7] employing a 
formation time of r = 1 fm/c. The latter quantity is uncertain by at least a factor of 
2 which implies a corresponding uncertainty in the energy density. Nevertheless, energy 
densities of a few GeV/fm^ suggest that the critical energy density for a QGP phase 
should be overcome in considerable space-time volumes at RHIC, where the relevant 
degrees of freedom are partons (quarks and gluons). Parton cascade calculations [8-10] 
are expected to provide suitable descriptions in the early phase [11,12] of these coUisions 
whereas hadrons should only be formed (by 'condensation') at a later stage which might 
be a couple of fm/c from the initial contact of the heavy ions. In fact, hybrid models like 
VNI+UrQMD [13], VNI+HSD [14] or the AMPT approach [15] also allow for a reasonable 
description of the 'soft' hadronic observables so far, which - due to the high interaction 
rate - are found to be close to the hydrodynamic limit [16]. On the other hand, once the 
local equilibrium limit is reached in the reaction, any conclusion on the dynamics in the 
early nonequilibrium phase and its dynamical degrees of freedom becomes highly model 
dependent. 

Moreover, the question of chiral symmetry restoration at high baryon density and/or 
high temperature is of fundamental interest, too [1,2]. Whereas lattice QCD calculations 
at zero baryon chemical potential indicate that a restoration of chiral symmetry goes along 
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with the deconfinement phase transition at a critical temperature Tc, the situation is less 
clear at finite baryon density where QCD sum rule studies show a linear decrease of the 
scalar quark condensate {qq) - which is nonvanishing in the vacuum due to a spontaneous 
breaking of chiral symmetry - with baryon density pb towards a chiral symmetric phase 
characterized by (qq) = 0. This decrease of the scalar condensate is expected to lead 
to a change of the hadron properties with density and temperature, i.e. in a chirally 
restored phase the hadrons might become approximately massless as suggested in Ref. 
[17]. However, chiral symmetry restoration only implies that vector and axial vector 
currents should become equal [18,19]. Thus vector and axial vector excitations of the 
QCD vacuum must have the same spectral functions in the chiral limit. As demonstrated 
in Refs. [20,21] such a restoration of chiral symmetry in central nucleus-nucleus collisions 
should — driven by the baryon density — occur at bombarding energies of 5-10 A-GeV. 
In Ref. [22] it has been argued, furthermore, that such 'phase transitions' should be seen in 
a much lower strangeness to entropy ratio. It has been also suggested [21] that especially 
the /7r~^ might give an indication for a chirally restored phase. 

The fact that the K~^/7r~^ ratio is found experimentally to be higher at top AGS 
energies of 11 A-GeV than at 160 A-GeV has raised speculations about the appearance 
of 'new physics' at energies between AGS and top SPS. To shed some light on this issue, 
the NA49 Collaboration has started an energy scan at the SPS. First results have become 
available now at 40 and 80 A-GeV [23-25] and further studies are foreseen at 30 and 20 
A-GeV [26]. Since this topic is of current interest we will restrict our investigations to the 
AGS and SPS energy range in this paper. 

Prom the theoretical side the various hadron spectra are conventionally calculated 
with nonequilibrium kinetic transport theory (cf. [27-32]). However, the calculated kaon 
to pion ratio from central nucleus-nucleus collisions turns out to vary by factors as large as 
2 if different transport approaches are applied [21,33,34,32]. Thus a unique interpretation 
of the data is questionable so far. On the other hand, statistical models [35] show a 
maximum of /-k^ ratio at ~ 30 A-GeV since the relative strangeness content of baryons 
is highest at low bombarding energies. It decreases with higher energies due to an increase 
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of temperature and a decrease of the baryon chemical potential. However, an analysis 
within the UrQMD transport model suggests that chemical and thermal equilibria are 
achieved only briefly in a small central overlap region of heavy-ion collision due to a very 
fast expansion of the hadronic fireball [36]. Moreover, the analysis of Ref. [37] within 
the HSD transport approach indicates that the equilibration time for strangeness at all 
bombarding energies is larger (> 40 fm/c) than the reaction time of nucleus-nucleus 
collisions. Thus the statistical model fits to the data have to be considered with some 
caution since they are not understood microscopically. 

In this work we concentrate on hadronic rapidity distributions of protons, kaons, an- 
tikaons and hyperons and their yields and ratios from Au + Au (or Pb + Ph) coUisions 
from SIS to SPS energies. The aim of our study is twofold: first, to find out the sys- 
tematic differences between two currently used transport approaches (denotes as UrQMD 
[38,39] and HSD [27,40]) and second, to look for common failures in comparison to re- 
lated experimental data that have become available recently [23-25] or provide predictions 
for experimental studies in the near future [26], which are also of relevance for the new 
GSI-proposal [41]. 

Our work is organized as follows: In Section 2 we will describe the main ingredients 
of the UrQMD and HSD transport approaches and point out conceptual differences. In 
Section 3 we study baryon stopping in central Au-\-Au collisions from 4 to 160 A-GeV in 
comparison to experimental data (whenever available) . Section 4 is devoted to a detailed 
comparison of both transport approaches on tt^, , and A+S° rapidity distributions, 
yields and different particle ratios as a function of bombarding energy from 2 to 160 
A-GeV. Again the calculations will be confronted with experimental data taken at the 
AGS and SPS. A direct comparison of UrQMD and HSD on the pp and 7r~p reaction 
level is given in Section 5 to quantify the differences in the 'elementary' differential cross 
sections. Section 6 concludes our study with a summary and discussion of open problems. 
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II. TRANSPORT MODELS - URQMD AND HSD 



In this work we employ two different transport models, i.e. the UrQMD and HSD 
approaches, that have been used to described nucleus-nucleus collisions from SIS to SPS 
energies for several years. Though different in the numerical realisation, both models are 
based on the same concepts: string, quark, diquark (g, q, qq, qq) and hadronic degrees of 
freedom. It is important to stress that both approaches do not include any explicit phase 
transition to a quark-gluon plasma (QGP). The philosophy is that a common failure of 
both models in comparison to experimental data should - model independently - indicate 
the appearance of 'new physics'. 

The UrQMD (Ultra-relativistic Quantum Molecular Dynamics) transport approach is 
described in Refs. [38,39]. It includes all baryonic resonances up to an invariant mass 
of 2 GeV as well as mesonic resonances up to 1.9 GeV as tabulated in the PDG [42]. 
For hadronic continuum excitations a string model (let's denote it as 'Frankfurt' string 
model (FSM)) is used. The hadron formation time (which relates to the time between 
the formation and fragmentation of the string in the individual hadron-hadron center-of- 
mass frame) is in the order of 1-2 fm/c depending on the momentum and energy of the 
created hadrons (using the "yo-yo" formation concept for the time definition) [38,39]. The 
UrQMD transport approach is matched to reproduce nucleon-nucleon, meson-nucleon and 
meson-meson cross section data in a wide kinematical regime [38,39]. At the high energies 
considered here the particles are essentially produced in primary high energy collisions by 
string excitation and decay, however, the secondary interactions among produced particles 
(e.g. pions, nucleons and excited baryonic and mesonic resonances) - that also contribute 
to the particle dynamics - are included as well. 

Whereas UrQMD operates as default in the cascade mode, i.e. with hadron potentials 
turned off, the HSD (Hadron- String Dynamics) transport approach includes (by default) 
scalar and vector fields of the particles which determine the mean-field propagation of 
the hadrons between collisions (cf. Fig. 2 of Ref. [21]). The HSD transport approach 
incorporates only the baryon octet and decuplet states and N*(1440), N*(1535) as well 
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as their antiparticles and the 0~ and 1~ meson octets. Higher baryonic resonances are 
discarded as exphcit states (for propagation) in HSD; they are supposed to "melt" in the 
nuclear medium even at normal nuclear density (see e.g. [43,44]). The argument here 
is that the resonance structure (above the A-peak) is not seen experimentally even in 
photoabsorption on light nuclei [45]. In contrast to the resonance concept - adopted in 
UrQMD for all low energy baryon-baryon and meson-baryon collisions - HSD includes 
the direct (non-resonant) meson production in order to describe the corresponding cross 
sections (for the details see Ref. [27]). 

In the HSD approach the high energy inelastic hadron-hadron collisions are described 
by the LUND string model (realized by FRITIOF-7.02 [46]), where two incoming hadrons 
emerge from the reaction as two excited color singlet states, i.e. 'strings' (as in UrQMD). 
The formation time of all hadrons — composed of light and strange quarks — in HSD is 
assumed to be rp ~ 0.8 fm/c in the hadron rest frame [27,40], which is lower than the 
'average' of the exponentially distributed formation times of 1-2 fm/c used in UrQMD. 
Note, that in both models the formation time in the calculational frame for heavy-ion 
collisions (laboratory or center-of-mass frame) is dilated by the Lorentz 7-factor, i.e. 

Since at high energy heavy-ion collisions particle production essentially proceeds via 
baryon-baryon and meson-baryon string excitations and decays, it is worth to discuss the 
differences in the reahzations of the string models used in UrQMD and HSD. In both 
string models the production probability P of massive ss or qqqq pairs is suppressed 
in comparison to light flavor production {uu, dd) according to a generalized Schwinger 
formula [47] 



with the string tension k ^ 1 GeV/fm. Thus in the string picture the production of 
strangeness and baryon-antibaryon pairs is controlled by the masses of the constituent 
quarks and diquarks. Inserting the constituent quark masses = 0.3 GeV and = 
0.5 GeV a value for the strangeness suppression factor 7^ ^ 0.3 is obtained. While 
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(1) 



the strangeness production in proton-proton collisions at SPS energies is reasonably well 
reproduced in the LUND string model with 7^ = 0.3, the strangeness yield for p + Be 
collisions at AGS energies (which is a good probe for the isospin averaged elementary p+p 
and p + n reactions) is underestimated by roughly 30% [33] . Therefore the strangeness 
suppression factor has been enhanced to 0.4 at AGS energies for the elementary nucleon- 
nucleon cross section in HSD. Thus, the relative production probabilities for the different 
quark flavours in the HSD model are 



HSD : u : d : s : diquark — < 



1 : 1 : 0.3 : 0.07 for > 20 GeV 

(2) 

1 : 1 : 0.4 : 0.07 for ^ < 5 GeV 



with a linear transition of the strangeness suppression factor as a function of s/s in 
between. The relative production probabilities for the different quark flavors in UrQMD 
are fitted to 

UrQMD : u : d : s : diquark = 1:1: 0.35 : 0.1. (3) 

Additionally fragmentation functions f{x,mt) must be specified, which are the probabil- 
ity distributions for hadrons with transverse mass rrit to acquire the energy-momentum 
fraction x of the fragmenting string. One of the most common fragmentation functions is 
used in the LUND model [46] (which is adopted in the HSD approach [33]) 

f{x,mt)^—{l — xyexp(—bmf/x], (4) 

with a = 0.23 and b = 0.34 GeV~^. In UrQMD different fragmentation functions are used 
for leading nucleons and newly produced particles, respectively (cf. Ref. [38] Fig. 3.16): 

/(a;)nuci = exp ^— ^ ^ , for leading nucleons 

/(x)prod = (1 ~ 2;)^, for produced particles (5) 

with A — 0.275 and B — 0.42. The fragmentation function /(x)prod — used for newly pro- 
duced particles — is the well-known Field- Fcynman fragmentation function [48]. At the 
string break-up the qq-pairs have zero transverse momenta in the string reference frame. 
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but the transverse momentum distributions of the single quark (pt) and the corresponding 
antiquark {—pt) are taken to be gaussian 



with a = 1.6 GeV/c. 

Despite the differences in the fragmentation functions, both string models describe 
quite well the data available for particle multiplicities and total spectra from pp collisions 
at high energies (see Ref. [38] for UrQMD and Ref. [33] for HSD). Also the inelastic pion- 
proton cross section is in good agreement with the experimental data in both models 
whereas differential spectra can differ substantially (cf. Section V). The LUND string 
model (in HSD) has also been tested for low energy pp collisions as well as for ttN inter- 
actions (cf. Chapter 2 in Ref. [33]). It has been shown that the (LUND) string model 
underestimates pion and kaon/antikaon yields closer to their production threshold. In 
HSD the threshold for string formation and decay thus is taken as ^/s = 2.6 GeV for 
baryon-baryon collisions and at y/s — 2.1 GeV for meson-baryon collisions. For lower in- 
variant energies a/s resonant and direct meson production mechanisms (e.g. ttN Ntttt) 
dominate, which are implemented in addition in HSD to ensure smooth excitation func- 
tions of the meson multiphcities from threshold to a few hundred GeV/c. 



Though various predictions have been made in both transport models since a couple of 
years, it is of importance to compare with actual data. We employ the experimental cuts 
in centrality to get as realistic as possible in the comparison of baryon stopping achieved 
in both theoretical approaches and in the different experiments. A related comparison is 
presented in Fig. 1 for protons from 5% (4, 6, 8, 10.7 and 160 A-GeV) and 7% central 



(20, 40, 80 A-GeV) Au + Au (AGS) and Pb + Ph (SPS) coUisions at 4-160 A-GeV ^ 




(6) 



III. BARYON STOPPING 



Note, that for all UrQMD and HSD calculations presented in this work the centrality of 
the reaction has been determined by a comparison of the transport calculations to the energy 
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The experimental data at 4, 6, 8, 10.7 A-GeV have been taken from Ref. [49] (circles), 
at 160 A-GeV from [50] (triangles) and from [51] (circles). The full symbols (here and 
for all further figures) correspond to the measured data whereas the open symbols are 
the data reflected at midrapidity. The solid hues with stars show the results from the 
UrQMD calculations while the sohd and dashed lines stem from the HSD approach with 
and without potentials, respectively. 

We note, that in the UrQMD calculations 'spectator' protons have been cut off whereas 
they are still present in the HSD calculations; this leads to the maxima in the proton 
rapidity distributions at target and projectile rapidity in the HSD calculations (cf. Fig. 1). 
Nevertheless, the HSD cascade calculations are found to agree with UrQMD cascade 
calculations from 4-20 A-GeV within 5%, whereas UrQMD shows somewhat more proton 
stopping than HSD at higher bombarding energies. The mean-field propagation effects 
in the HSD approach are most pronounced at low bombarding energies leading to a 
reduction of baryon stopping from 4 - 10.7 A-GeV and a flatter rapidity distribution 
dN/ dy around midrapidity slightly closer to the experimental data. This effect can easily 
be attributed to the energy stored in the repulsive mean field at high baryon density and 
moderate bombarding energy. Above about 40 A-GeV such potential effects are no longer 
statistically significant in the calculations since the repulsive mean field decreases strongly 
with momentum (cf. Fig. 2 of Ref. [21]) such that at 40 A-GeV no repulsion is seen by the 
nucleons in the initial high density phase. Only when the system partly thermalizes and 
the nucleon momenta relative to the fireball reference frame become smaller the baryons 
'feel' again a repulsive mean field, however, now at rather low baryon density. We recall 
that at density po the potential is even attractive for momenta < 600 MeV/c. 

In general the HSD results indicate slightly less baryon stopping than the UrQMD 
calculations. At 160 A-GeV the experimental data favor a minimum of the distribution 



distribution in the Veto-calorimeter of the NA49 collaboration for SPS energies and to the New 
Multiplicity Array (NMA) and ZCAL-calorimeters at AGS energies 



10 



at midrapidity, which is reproduced by the HSD calculations. However, the UrQMD 
calculations only deviate by ~ 5%. We note, that for semi-central and peripheral Pb + Pb 
collisions at 160 A-GeV the UrQMD calculations are in good agreement with the data 
from Ref. [50] (cf. Fig. 1 in Ref. [52]). 

Thus, the overall description of the rapidity spectra from both models (with/without 
potentials) in this wide energy regime is quite remarkable in view of the different 'hadronic' 
degrees of freedom and string 'parameters' involved. 

IV. PION AND STRANGENESS PRODUCTION 

We continue with tt^, K^, K~ and A + E° rapidity spectra in 5%, 7%, or 10% central 
collisions of Au-\-Au or Pb+Pb, respectively, from 4-160 A-GeV. We compare to the AGS 
data from Refs. [53-58] in Fig. 2 and SPS data from Refs. [23-25] in Fig. 3. Here the 
thick solid lines denote the HSD results including the potentials, the dashed lines represent 
HSD calculations in the cascade mode, which should be directly compared to the UrQMD 
results (thin solid lines with stars). Both Figs. - taken together - provide an overview 
on the energy dependence of the different rapidity distributions and the virtues/failures 
of the transport models. 

At 4 A-GeV all transport versions overestimate the 7r+ spectra. This deviation is most 
pronounced in the HSD 'cascade' version. On the other hand, UrQMD is higher than 
HSD in the strangeness channels K"^ and A-|-E°, whereas HSD (with potentials) is quite 
compatible with the midrapidity data. At 6 A-GeV the HSD (with potential) calculations 
are in line with the experimental K"^ and A + E° data whereas the UrQMD results are 
still too high. All calculations again overestimate the tt"*" rapidity spectra. At 8 A-GeV 
the UrQMD and HSD results continue to overestimate the tt"*" yield, while now the HSD 
calculations (with potential) fall low in the strangeness channels contrary to the HSD 
cascade and UrQMD results. At 10.7 A-GeV the picture continues: both models give too 
many tt"^, the strangeness yield is approximately reproduced in the cascade versions and 
underestimated in the HSD simulation with potentials. 
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At 20 A-GeV, most relevant for the future GSI facility [41], both transport models 
(with/ without potentials) give the same rapidity distributions for all hadrons considered 
here. Note, that the cascade HSD results are not shown explicitly in Fig. 3 since the 
cascade and potential results are practically coincident. As discussed above this is due to 
the momentum-dependence of the scalar and vector potentials in HSD; both potentials 
vanish for relative momenta of a couple of GeV/c. Deviations between HSD and UrQMD 
come up again at 40 A-GeV, where UrQMD now gives more n' and lacks a few K'^, 
whereas the K" and A + E° spectra are reasonably reproduced. At this energy HSD 
(with/without potentials) only overestimates the 7r~ distribution slightly. At 80 A-GeV 
the UrQMD spectra are also fine for K~ and A, but low for K'^ and too high for tt". At 160 
A-GeV the HSD calculations indicate slightly too few 7r~ at midrapidity, however, do well 
for K"^ and A + S° as already demonstrated earlier in Refs. [27,33] . The UrQMD results 
here show too many tt"^ and a shghtly higher A + E'^ yield, whereas the K'^ distribution 
is underestimated. We note, that the UrQMD results presented in Fig. 3 are taken from 
Ref. [34]. 

As a more general overview on the tt^ abundancies in central Au + Au and Ph + Ph 
collisions we show in Fig. 4 the tt"*" (upper plots) and tt" (lower plots) multiplicities at 
midrapidity (left column) and integrated over rapidity (right column) as a function of the 
bombarding energy in comparison to the available data from Refs. [23,53]. Here the solid 
hues with open triangles show the results from the UrQMD calculations while the solid 
lines with open squares and dashed lines stem from the HSD approach with and without 
potentials, respectively. At lower AGS energies the UrQMD model gives slightly less pions 
then HSD (with/without potential), but both models overpredict the midrapidity data 
(except UrQMD at 2 A-GeV, which is in line with data point at midrapidity). About 20 
A-GeV is the "crossing point" for both transport calculations and at SPS energies the 
tendency turns around: UrQMD gives more pions than HSD, so that HSD is now in a 
better agreement with the experimental data. 

The overestimation of the pion yield at low and high energies by the transport models 
is presently not well understood. In Ref. [59] Larionov et al. have speculated that higher 
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baryon resonances - more massive than the A-resonance - might be quenched at densities 
above ~ 1.5 po both in A^'A'" as well as ttA^ channels. The net effect at SIS energies (and 
slightly higher bombarding energies) is a reduction of the pion yield in collisions of heavy 
systems such as Au+Au whereas light systems (e.g. C+C) are not effected very much [59]. 
Though this might be a possible explanation at SIS and lower AGS energies, other effects 
such as strong pion self energies could also lead to a reduction of the pion abundance in the 
transport models. Such many-body effects have not been incorporated in the calculations 
presented here. At SPS energies and especially at 160 A-GeV another production channel 
for pions becomes sizeable, i.e. the annihilation of baryon-antibaryon {BB) pairs that 
leads on average to 5 pions (or more). However, by employing detailed balance on the 
many-body level [60] the BB annihilation rate is found to be almost compensated by 
many-meson fusion channels that recreate BB pairs. This many-body channel is not 
incorporated in UrQMD and thus might partly be responsible for the overestimation of 
the pion yield. A rough estimate, however, shows that this missing channel cannot be the 
only reason: as calculated in Ref. [60] (and approved by recent NA49 measurements [25]) 
thep/ < TT > ratio (< tt >= 0.5(77"*" -|-7r~)) is approximately independent on the centrality 
of the Au + Au collision and ~ 1%. Assuming the number of n to be equal to the number 
of antiprotons as well as the number of antihyperons [25] and adopting isospin symmetry 
for the pions, we get N/n 1 %. If all TV's produced initially are annihilated, which 
is an upper estimate for the pion production from this channel, then we could achieve a 
maximum increase of the pion number by 5 %. Thus a ~20 % overestimation of the pion 
yield by UrQMD at 160 A-GeV can only partly be attributed to the missing multi-meson 
fusion channels. 

The (upper plots) and K~ (lower plots) multiphcities at midrapidity (left column) 
and integrated over rapidity (right column) are shown in Fig. 5 as a function of the 
bombarding energy in comparison to the available data from Refs. [23,53]. Here the 
midrapidity and total yields summarize the findings from Figs. 2-3: The abundancies 
are well described by both transport models. The yield at midrapidity (left column) is 
slightly overestimated by UrQMD at AGS energies and underestimated at SPS energies, 
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whereas HSD is in a reasonable agreement with the data except of the upper AGS energies 
(with potentials included). This tendency stays the same for the Att kaon yield, however, 
at SPS energies the difference between both models for the iir yields is smaller than 
at midrapidity since UrQMD provides a slightly broader kaon rapidity distribution than 
HSD (cf. Fig. 3). Thus, an underestimation of strangeness production is not the prevaihng 
issue as demonstrated in Fig. 5 in comparison to the recent experimental data from NA49 
[23]. Both transport models can roughly describe - within their systematic range of 
uncertainties - the spectra and abundancies. 

In Fig. 6 we show the /t:^ and K~ /t{~ ratios at midrapidity (left column) and 
integrated over all angles (right column) as a function of the bombarding energy for central 
collisions of Au + Au (AGS) or Ph + Ph (SPS) in comparison to the available data from 
Refs. [23,53]. Whereas the excitation function of the K~ /ti' ratio is roughly reproduced 
by both transport models, the maximum in the j'K^ ratio seen experimentally both at 
midrapidity (upper left part) and in 47r (upper right part) is not described by HSD as well 
as UrQMD. For the /■n'^ ratio both models give quite different results. HSD gives a 
monotonous increase of this ratio with bombarding energy (as pointed out in Refs. [21,33] 
^), whereas within UrQMD the ratio shows a maximum around 10 A-GeV and then drops 
slightly for the midrapidity ratio or stays roughly constant for the 47r ratios. In view of 
Figs. 2-5 this failure is not primarily due to a mismatch of strangeness production, but 
more due to an insufficient description of the pion abundancies. 

Fig. 7 shows the excitation functions of A + S°-hyperons at midrapidity (upper left 
part) and integrated over all angles (upper right part) as a function of the bombarding 
energy for central collisions of Au + Au (AGS) or Ph + Ph (SPS) in comparison to the 
available data from Refs. [24,25,57,58,61]. Here UrQMD (solid lines with open triangles) 

^We note that the HSD results presented in this work are produced with much higher statistics 
than in Refs. [21,33] due to the increasing computer power available. Also the centrality selection 
is done now in line with the actual experimental set-up. 
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slightly underestimates the 47r yields at 40 and 80 A-GeV whereas HSD (solid hnes with 
open squares) seems to give a better description. Nevertheless, all models compare rather 
well with data. 

The (A + E°)/7r ratios^ at midrapidity (lower left part) and integrated over An (lower 
right part) are underestimated slightly which should again be attributed to the pion 
excess in the transport models (see above). Nevertheless, the maxima in the ratios (47r 
and midrapidity) observed experimentally is qualitatively reproduced by both models 
indicating that with increasing bombarding energy s-quarks are more frequently produced 
within mesons {K, K*) rather than in associate production with baryons. A similar trend 
is also found in statistical model fits [35]. 

The excitation function for the K~ / K'^ ratio in central collisions of Au + Au (AGS) 
or Ph + Ph (SPS) is shown in Fig. 8 for midrapidity ratios (l.h.s.) and in 47r (r.h.s.). 
Experimental data [23,53] here are only available for the midrapidity ratios. Again we 
find that both transport models give similar results for this ratio which are comparable 
to the data. Statistical models also fit this ratio quite well. Whether this finding implies 
that 'chemical equilibration' is approximately achieved in central collisions of heavy nuclei 
is still an open question [62]. 

One note of caution again: the pions in both transport models are treated as 'free' 
particles, i.e. with their vacuum mass. On the other hand, lattice QCD calculations as 
well as effective Lagrangian approaches like the Nambu-Jona-Lasinio (NJL) model show 
that the pion mass should increase with temperature and density. So the overestimation 
of the pion yields in HSD and UrQMD might be a signature for a dynamically larger 
pion mass. Moreover, in-medium changes of the strange hadron properties, as known 
from experimental studies at SIS energies, should also show up at AGS and SPS energies. 
Thus, including all medium effects simultaneously in a consistent way might provide a 
more conclusive interpretation of the ratios in Figs. 6-8. However, such calculations require 



^The pion multiplicity is calculated as tt = 3/2(7r+ -|- tt ) in line with Ref. [24,25] 
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a precise knowledge about the momentum and density dependence of the hadron self- 
energies which is not available so far. Note, that up to now in-medium modifications of the 
iiT"*", K~ properties have been studied with HSD employing a dropping of K~ and increase 
of K'^ masses in the medium [63]. As summarized in [27] such a scenario leads to an 
enhancement of K~ and a lowering of yields at SIS energies (which is close to threshold 
for K~^,K~ production). It modifies only slightly the strangeness abundancies at SPS 
energies. However, chiral symmetry restoration also requires a simultaneous modification 
of the pion properties. 

We close this Section with some speculations about the failure of the transport mod- 
els for the K~^/7r~^ ratio at the top AGS energies (and slightly above). To this aim let's 
assume that the pion yield is decreased by some mechanism to the actual yield observed 
experimentally. Since kaons and antikaons are also produced in secondary non-strange 
meson-baryon collisions, this will imply also a reduction of the kaon number in the trans- 
port calculations. For central Au+Au reactions these secondary channels give roughly the 
same amount of kaons and antikaons as the primary NN channels. Thus a 20% decrease 
in the pion number will lead to a maximum change in the kaon number by 10% and the 
K'^/n'^ ratio might increase by about 10%, too. This relative increase will improve the 
situation in comparison to experiment, however, not solve the problem. Another possibil- 
ity is the enhanced production of ss pairs in the hot and dense hadronic medium relative 
to the vacuum, i.e. an increase of 7^ in Eq. (1). Such modifications might be driven by 
an enhanced string tension k in the medium [64] or a relative decrease of mj. — rrig with 
density and temperature. Alternatively, also the formation of a QGP in the initial phase 
might lead to enhanced strangeness production. According to the authors point of view 
such phenomena cannot be excluded presently, however, there is also no strong evidence 
in favor of them. 
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V. URQMD VERSUS HSD 



Though both transport approaches - HSD and UrQMD - give quahtatively and quan- 
titatively similar results for proton, tt^, and hyperon rapidity distributions, there are 
quantitative differences that should be cleared up. To this aim we concentrate on idealized 
model comparisons which are not directly comparable to experimental data. 

The excitation functions for s- and s-quark production in 6 = collisions of Au + Au 
(AGS) or Pb + Pb (SPS) nuclei (HSD: l.h.s.; UrQMD: r.h.s.) are shown in Fig. 9 which 
finally end up in mesons or baryons (antibaryons). The solid (dashed) lines with full (open) 
squares and circles indicate the number of s-quarks (s-quarks) in baryons (antibaryons) 
and mesons, respectively. Both models give qualitatively and quantitatively similar results 
(as already found from the previous comparisons) showing that at ~ 80 A-GeV the same 
amount of s-quarks end up in mesons and baryons, whereas baryons are preferred at 
lower energies. The s-quarks in antibaryons are almost negligible in this energy range. 
We mention here that strangeness conservation is exactly fulfilled in both transport models 
such that at all energies the number of s-quarks is identical to the number of s-quarks. 

The channel decomposition (fraction in %) for the finally observed (left column) 
and K" (right column) from HSD (upper part) and UrQMD (lower part) are shown in 
Fig. 10 as a function of bombarding energy for Au + Au (AGS) or Pb+Pb (SPS) coUisions 
at impact parameter 6 = 0. Note, that only the contributions of the dominant channels 
are shown in Fig. 10. For a correct interpretation of Fig. 10 one has to keep in mind, that 
at high energies initially s, s quarks arc produced in primary nucleon-nucleon collisions 
and later on in meson-baryon interactions via string excitation and decay. However, 
subsequently the strange particles (produced initially) participate in chemical reactions 
with flavor exchanges. As a result only a few percent of the 'primary' kaons/antikaons 
remain unaffected by secondary inelastic interactions (cf. the lines denoted as ^BB string' 
in Fig. 10). Most of the final K'^ and K~ mesons - above 10 GeV - finally stem from 
K*^{892) decays (fines 'i^*(892) decay') which are either produced directly in string 
decays or by pion-kaon resonant scattering. As seen from the lower part of Fig. 10, 
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in UrQMD a large fraction of final K'^ ^ K~ stem directly (without K* production and 
decay) from meson-baryon collisions (line 'mS string'). These are realized in UrQMD via 
an excitation of a single string that furtheron decays isotropically, which is reminiscent of 
the resonance mechanism. The same mechanism is also used for high energy meson-meson 
collisions (hne 'mm string'). In HSD this channel is treated differently, i.e. via ir*(892) 
resonance production, and thus contributes to the 'i^*(892)' channel. Note also, that in 
UrQMD - in the channel denoted as 'mS string' - the kaon/antikaon-baryon coUisions 
are counted, too, whereas they are not counted here in HSD. In both models only a few 
percent of the final and K~ appear from the 0(1020) meson decays (lines '0 decay'). 
We mention that only a small fraction of the 0-decays can be reconstructed from K'^K~ 
invariant mass spectra. 

The conceptual differences in the treatment of strangeness production in both trans- 
port approaches are more pronounced at low energies. UrQMD implements the full reso- 
nance dynamics assuming the vacuum resonance properties: At low energies the baryon 
and meson resonances are excited in baryon-baryon (BB), meson-baryon {mB) or meson- 
meson (mm) coUisions (below the string thresholds) and decay (after time l/F) to K^, K' 
(cf. the line 'A?"*' in Fig. 10 denoting the fraction of the final K'^ ^ K~ from strange baryon 
resonance decays: iV*(1650), iV*(1710), iV*(1720), iV*(1990) and line 'F*' for the higher 
hyperon resonances as well as the lines 'Xi -|- and ao + /o indicating the fraction from 
the decay of meson resonances: iri(1270), Kq{1A?,{)) and ao(980), /o(980), respectively). 
In HSD these high strange baryon and meson resonances are not produced and prop- 
agated explicitly (as indicated in Section 2). Instead - in low energy BB^mB or mm 
collisions - strangeness is directly produced with respect to the corresponding transition 
rates and 2, 3 or 4-body kinematics for the final states (cf. lines denoted as 'vriV — >■ KY\ 
'ttF — > KN^ or 'mm coll.' in the upper part of Fig. 10). This comparison demonstrates 
that the individual channels are treated quite differently though the final results are very 
similar. 

Thus, as seen from Fig. 10, only a small fraction of kaons/antikaons (less than 10% 
of the kaons and less than 6% of the antikaons) from energetic initial baryon-baryon 
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collisions (cf. lines denoted as BB string) survives the hadronic rescattering phase during 
the expansion of the fireball without reinter action. Most of the final strange particles 
emerge after rescattering — shifting s quarks from mesons to baryons and vice versa — 
thus providing a very distorted picture of the initial strangeness production mechanism 
and the elementary degrees of freedom involved. Consequently, as pointed out in Ref . [34] , 
the K"^ and A spectra do not allow for stringent conclusions on the initial phase of high 
energy density. On the other hand, these frequent flavor exchange reactions may be 
viewed as the reason why statistical models (employing chemical equilibrium) seem to 
work reasonably well. 

A. pp reactions 

In order to get some idea about the differences between both transport approaches we 
go back to the description of the elementary channels like pp or 7T~p in vacuum. In this 
respect we show in Fig. 11 the proton rapidity distributions for pp collisions from HSD 
(solid lines) and UrQMD (dashed lines) between 4 and 160 GeV laboratory energy. This 
provides information on the different string excitation and fragmentation schemes. The 
experimental data at 160 GeV are taken from Ref. [50]. As seen from Fig. 11, however, 
the differences between the two string fragmentation schemes are only minor. Where do 
the differences in baryon stopping — shown in Fig. 1 — come from? 

The differential rapidity distributions for n^^K^ and A + E°'s from pp collisions, 
however, show substantial differences as demonstrated in Fig. 12. The experimental data 
for i^+'s and A + E°'s at 160 GeV are taken from Refs. [65] and [66], respectively. Though 
the TT rapidity distributions are roughly comparable in both models, there is a trend 
for UrQMD to give slightly more pions with decreasing bombarding energy than HSD, 
whereas in heavy-ion collisions the trend is opposite — HSD gives more pions at low 
energies than UrQMD, whereas UrQMD gives more pions at high energies^ (cf. Figs. 2, 



^We will attribute this to the influence of secondary mB scatterings, where more mesons are 
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3). For K~ mesons the results of both models are comparable at high energies, however, 
they deviate closer to the K~ production threshold. For and A+E° both models differ 
substantially, too. Here HSD gives more at low energies whereas the UrQMD rapidity 
distribution is broader at 160 GeV. The A yield from pp coUisions is also higher from HSD 
- due to strangeness conservation - and shows distinct peaks in the rapidity distribution 
closer to target and projectile rapidities at high energies, whereas the UrQMD rapidity 
distributions for A's are narrower and almost peaked at midrapidity. Experimental data 
[66] - available at 160 GeV - show a minimum at midrapidity giving no strict preference 
for one of the string fragmentation schemes. 

Thus strange quarks are produced more at midrapidity in UrQMD, both for mesons 
and baryons, whereas in HSD s-quarks are concentrated in mesons at midrapidity and 
in baryons at larger rapidities. These differences in the elementary differential rapidity 
spectra explain also the different rapidity distributions from central nucleus-nucleus col- 
hsions in Figs. 2 and 3 to a large extent. Presently it is not clear - due to the lack of 
corresponding experimental data - which fragmentation scheme is 'more realistic'. On the 
other hand, this comparison sheds some light on the 'systematic' uncertainties in present 
relativistic transport approaches. These 'systematic' uncertainties have to be kept in 
mind when attempting to draw conclusions from nucleus-nucleus collisions in comparison 
to experimental data. 

Before closing this Subsection we confront both transport models with the available 
data on the production cross sections of pions and strange hadrons from pp collisions. 
In Fig. 13 we show the inclusive tt"*", 7r~, , K~ and A -|- E° production cross sections 
from pp collisions versus the kinetic energy Eiah- The solid lines with open triangles 
show the UrQMD results, the solid line with full squares indicate the HSD results with 
the strangeness suppression factor 7^ defined by Eq. (2), whereas the dot-dashed lines 
correspond to 7^ = 0.3. The experimental data are taken from Refs. [67] (full triangles). 



produced within UrQMD. 
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[68] (full and open circles), [69] (open squares) and [70] (stars). The pion cross sections 
are quite well described by both models; UrQMD gives more 7r+ than HSD at low energy 
in line with the data point from [67], whereas HSD follows more closely the data from 
Ref . [68] . The inclusive antikaon cross section is well reproduced by both approaches. As 
already demonstrated in Fig. 12 HSD gives more K'^ and neutral strange hyperons than 
UrQMD below Eiab ~ 80 GeV. The neutral hyperon yield from UrQMD (for Eiab < 80 
GeV) is more in line with the data, whereas the K'^ yield is slightly underestimated from 
10 - 80 GeV. In contrast HSD seems to better reproduce the K'^ cross sections but to 
overestimate the A + E° yields in pp reactions at lower energies. 

The differences in these 'in-put' cross sections are quite sizeable, however, one has to 
keep in mind that only a single isospin channel is probed in Fig. 13 in comparison to 
data, whereas in nucleus-nucleus collisions essentially isospin averaged cross sections are of 
relevance. In fact, both transport models differ in the isospin dependent cross sections for 
A'^A^ collisions, whereas isospin averaged particle yields are more similar. We recall again 
that strangeness conservation holds explicitly for both transport models with respect to 
all reactions employed. 



B. vr p reactions 

We now turn to the elementary pion-nucleon collisions that play a substantial role in 
secondary meson-baryon coUisions. The differential tt"^, and A(-|-E°) rapidity distribu- 
tions from 7r~p reactions from 2-8 GeV/c laboratory momentum are shown in Fig. 14 for 
UrQMD (dashed lines) and HSD (solid lines). Also here experimental data are not avail- 
able for a comparison. Though the total and elastic Tr~p cross sections are very similar in 
both models and in line with experimental data (cf. [33,38]), the explicit rapidity distri- 
butions for various final states differ by up to a factor of 2-3. This holds true also for the 
isospin dependent cross sections (e.g. vs. K^) that are not probed in nucleus-nucleus 
collisions due to initial isospin averaging. In general the string model in UrQMD produces 
substantially more 7r~^,K^ etc. in tt -|- A^ reactions than the LUND model employed in 
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HSD. Consequently, hadron (including strange meson) production by secondary meson- 
baryon collisions is sizeably higher in UrQMD than in HSD. This observation clarifies to 
some extent the higher tt^ yield in Fig. 4 from UrQMD at SPS energies relative to HSD 
and the experimental data. On the other hand, strangeness production A + E°) from 
pp coUisions is much higher in the LUND approach (cf. Fig. 12) than in the FSM (used in 
UrQMD) such that one might expect the HSD approach to give more kaons and hyperons 
in central nucleus-nucleus collisions due to a higher initial production. As seen from Figs. 
5-7 this expectation does not hold true since in UrQMD the strangeness production from 
secondary {mB) channels is substantially higher now, which compensates - relative to 
HSD - for the initially lower strangeness production from NN collisions. 

In summary, the dominant differences between HSD and UrQMD for central nucleus- 
nucleus collisions can be traced back to difi^erent string fragmentation schemes for BB 
and mB strings that lead to substantially different hadron distributions in rapidity as 
well as isospin. Presently, these string models are not sufficiently controlled by differential 
experimental data. Furthermore, the string models employed are not tailored to describe 
the isospin dependence of the elementary cross sections at lower invariant energies ^/s. 

C. pp versus central AA reactions 

In order to explore the main physics from central AA reactions it is instructive to 
have a look at the various particle multiplicities relative to pp collisions as a function of 
bombarding energy. To this aim we show in Fig. 15 the total multiphcities of tt"*", K'^ and 
K~ (i.e. the 47r yields) from central Au + Au (at AGS) or Ph + Pb (at SPS) collisions 
in comparison to the total multiplicities from pp collisions versus the kinetic energy per 
particle Eiab- The solid fines with full triangles and squares show the UrQMD (l.h.s.) 
and HSD results (r.h.s.) for AA collisions, respectively, while the dotted lines with open 
triangles and squares correspond to the pp multiphcities calculated within UrQMD (l.h.s.) 
and HSD (r.h.s.). The multiplicities frompp reactions in Fig. 15 have been multiplied by a 
factor of 350/2 which corresponds to the average number of participants Apart in the heavy- 
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ion reactions for the centrality class considered divided by the number of participants in 
the pp reaction. We mention, that the comparison at lower bombarding energies of 2-4 
A-GeV has to be taken with some care due to the different influence of Fermi motion - in 
case of AA reactions - on the production of pions and K"^ mesons. 

The general trend from both transport approaches is quite similar: we observe a 
slight absorption of pions at lower bombarding energy and a relative enhancement of 
pion production in heavy-ion collisions above about 10 A-GeV. Kaons and antikaons from 
A A collisions are always enhanced in central reactions relative to scaled pp multiplicities. 
This enhancement is more pronounced within UrQMD than in HSD due to the larger cross 
sections employed in ttN secondary reactions as demonstrated in the previous Subsection. 

VI. SUMMARY 

In this work we have performed a systematic analysis of hadron production in central 
Au + Au or Ph + Ph collisions from SIS to SPS energies within the HSD (with and 
without potentials) and UrQMD transport approaches in comparison to the experimental 
data available. We flnd that both transport approaches — which are based on quite 
different initial ingredients — roughly give comparable results for the different p, tt^, 
and hyperon distributions in a wide energy regime from 2 - 160 A-GeV. It is remarkable 
that the cascade mode of HSD (which operates by default in the potential mode) gives 
to a large extent comparable results for strangeness production as the UrQMD cascade. 
This observation suggests that — inspite of the different elementary 'input' cross sections 
— the systems might reach approximate chemical equilibrium. This is a prerequisite for 
an analysis within statistical models [35,62,71]. In fact, the channel decomposition of 
strangeness production chains in both models are quite different (cf. Fig. 10) since the 
degrees of freedom (hadron resonances and strings) substantially differ for collisions at 
hadron-hadron coUision energies around 2-3 GeV in the region of string thresholds. 

We have found that at SPS energies HSD and UrQMD quite well reproduce the ex- 
perimental data for K~ and A -|- E° rapidity distributions at midrapidity as well as the 
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An yields. At 20 A-GeV both models agree very well among each other for all hadrons. 
This provides rather solid predictions for the future GSI heavy-ion program [41]. At 
AGS energies (< 11 A-GeV) the yield is slightly overestimated by UrQMD (except 
of 10.7 A-GeV), whereas HSD underestimats kaon production at the upper AGS energies 
(especially with baryon potentials included) . The K~ and A + E° data are reasonably 
described by both models. We have found also that HSD and UrQMD differ in the pion 
multiphcities — at lower AGS energies the UrQMD model gives slightly less pions than 
HSD (with/without potential), but both models overpredict the midrapidity data (except 
UrQMD at 2 A-GeV). At SPS energies the tendency turns around: UrQMD gives more 
pions than HSD, such that HSD is now in a better agreement with the experimental data. 
These differences between the transport approaches could be traced back to a large extent 
to different string fragmentation schemes which presently are insufficiently controlled by 
experimental data at the energies of interest here. 

The excitation functions of pions, kaons and antikaons from central Au + Au (or 
Ph + Pb) collisions relative to scaled pp reactions from the two transport models are very 
similar: both approaches give an absorption of pions at lower bombarding energy and a 
relative increase of pion production for Eiab > 10 A-GeV. Kaons and antikaons from AA 
collisions are enhanced in central reactions relative to scaled pp collisions at all energies 
by a factor of > 2. 

We have found that the failure of both models to reproduce the experimental excitation 
function for the /ti^ ratio in central nucleus-nucleus collisions — which might suggest 
the presence of a different state of hadronic matter in the early phase of these collisions 
— is not primarily due to an underestimation of strangeness production. Our systematic 
study in comparison to the most recent data from the NA49 Collaboration demonstrates 
that this failure is mainly due to an inadequate description of pion dynamics. We attribute 
this to the fact that the pions in both transport models are treated as 'free' on-shell 
particles, i.e. with their vacuum properties and 5-like spectral functions in mass. On the 
other hand, lattice QCD as well as effective Lagrangian models indicate an increase of the 
pion mass with temperature and density. Furthermore, the pion spectral function should 
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become broad in the medium due to the interactions. All these medium modifications 
have not been included in the calculations presented in the work. Thus the overestimation 
of the pion yields could be a signature for a chiral symmetry restoration which might occur 
at the high baryon/meson densities achieved in relativistic heavy- ion collisions. Including 
the medium effects for pions and all strange particles simultaneously in a consistent way 
in an 'off-shell transport approach' [72] could provide a more conclusive interpretation of 
the experimental data. This, however, requires a precise knowledge about the momentum 
and density dependence of the hadron self-energies in a wide energy regime and full off- 
shell transition matrix elements [72] . Such a program is clearly beyond the scope of our 
present study. 

Another problem of the transport approaches used here is that detailed balance is not 
implemented for n m transitions with n,m > 2 [36]. Thus multi-particle collisions 
might change the dynamical picture accordingly and lead to 'shorter' chemical equilibra- 
tion times [60,73,74]. In fact, the importance of 3 2 transitions has been demonstrated 
in the extended HSD transport approach in Ref. [60] for antibaryon reproduction by me- 
son fusion for 74-1-^4 coUisions at the AGS and SPS. In order to achieve a more conclusive 
answer from transport studies multiparticle interactions will have to be included in future 
generations of transport codes. 

What to conclude from the detailed comparisons presented in this work? Coming 
back to the question raised in the introduction about common failures in comparison to 
related experimental data, we can quote an insufficient accuracy in the description of 
the pion degrees of freedom by both transport models. Does this provide a signal for 
'new physics' in view of a QGP? The answer of the authors to this question with respect 
to the experimental observables studied is: no! As discussed above, the 'systematic 
uncertainties' in the 'on-shell' transport approaches are within the range of the deviations 
seen in comparison to the data or even larger. Furthermore, the question raised in the 
title of this paper - anything strange with strangeness? - also has to be answered with 
'most likely not' ! 
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y y 

FIG. 1. The rapidity distributions for protons from 5% (4, 6, 8, 10.7 and 160 A-GeV) and 7% 
central (20, 40, 80 A-GeV) Au + Au (AGS) and Pb + Pb (SPS) collisions at 4-160 A-GeV. The 
experimental data at 4, 6, 8, 10.7 A-GeV have been taken from Ref. [49] (circles), at 160 A-GeV 
from [50] (triangles) and from [51] (circles). The full symbols correspond to the measured data, 
whereas the open symbols are the data reflected at midrapidity. The solid lines with stars show 
the results from the UrQMD calculations while the solid and dashed lines stem from the HSD 
approach with and without potentials, respectively. 
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FIG. 2. The rapidity distributions for tt^jK^, K~ and A + S^'s from 5% central Au + Au 
collisions at 4-10.7 A-GeV in comparison to the experimental data from Refs. [53-58]. The thin 
solid lines with stars show the results from the UrQMD calculations while the thick solid and 
dashed lines stem from the HSD approach with and without potentials, respectively. 
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FIG. 3. The rapidity distributions for tt~,K+,K~ and A + S°'s from 5% (160 A-GeV), 7% 
(20, 40 and 80 A-GeV) or 10% central (A + S° at 160 A-GeV) Pb + Pb collisions at 20-160 
A-GeV in comparison to the experimental data from Refs. [23-25]. The solid lines with stars 
show the results from the UrQMD calculations while the thick solid and dashed lines stem from 
the HSD approach with and without potentials, respectively. 



34 



250 



200 



150 



50 


250 

200 

O 

^ 150 

Z, 100 

50 



71 




71 



# exp. data 
-□ — HSD, pot. 

HSD, cas. 

-A— UrQMD 



10 

"T 



100 

n-ri 




# exp. data 
-□ — HSD, pot. 

HSD, cas. 

-A— UrQMD 



10 

EJA [GeV] 



100 



800 



600 
^ 400 



200 


800 



i 1 ll 




' ' "1 1 ' ' ' 

A 




+ 


/ 




71 


A 


- 




/A 

d/ 












y 9 exp. data 






^ — □ — HSD, pot. 






^ HSD, cas. 


. . .1 




—A— UrQMD 






# exp. data 
—□—HSD, pot. 

HSD, cas. 

—A— UrQMD 



10 100 

E,7A [GeV] 



FIG. 4. The excitation function of vr"'" (upper plots) and vr (lower plots) yields from 5% 
(AGS energies and 160 A-GeV) or 7% central (20, 40 and 80 A-GeV) Au^Au (AGS) or Pb + Pb 
(SPS) collisions in comparison to the experimental data from Refs. [53,23] for midrapidity (left 
column) and rapidity integrated yields (right column). The solid lines with open triangles show 
the results from the UrQMD calculations while the solid lines with open squares and dashed 
lines stem from the HSD approach with and without potentials, respectively. 
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FIG. 5. The excitation function of (upper plots) and K (lower plots) yields from 5% 
(AGS energies and 160 A-GeV) or 7% central (20, 40 and 80 A-GeV) Au + Au (AGS) or Pb + Pb 
(SPS) collisions in comparison to the experimental data from Refs. [53,23] for midrapidity (left 
column) and rapidity integrated yields (right column). The solid lines with open triangles show 
the results from the UrQMD calculations while the solid lines with open squares and dashed 
lines stem from the HSD approach with and without potentials, respectively. 
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FIG. 6. The excitation function of the K~^/tt~^ (upper plots) and K /it (lower plots) 
ratios from 5% (AGS energies and 160 A-GeV) or 7% central (20, 40 and 80 A-GeV) Au + Au 
(AGS) or Pb + Pb (SPS) collisions in comparison to the experimental data from Refs. [53,23] 
for midrapidity (left column) and rapidity integrated yields (right column). The solid lines with 
open triangles show the results from the UrQMD calculations while the solid lines with open 
squares and dashed lines stem from the HSD approach with and without potentials, respectively. 
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FIG. 7. The excitation function of the A + S*^ yields and (A + S'^)/7r ratio from 5% (AGS 
energies), 7% (20, 40 and 80 A-GeV) or 10% central (160 A-GeV) Au + Au (AGS) or Pb + Pb 
(SPS) collisions in comparison to the experimental data from Refs. [24,25,57,61] for midrapidity 
(left column) and rapidity integrated yields (right column). The solid lines with open triangles 
show the results from the UrQMD calculations while the solid lines with open squares and 
dashed lines stem from the HSD approach with and without potentials, respectively. 
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FIG. 8. The excitation function of the K /K^ ratio from 5% (AGS energies and 160 A-GeV) 
or 7% central (20, 40 and 80 A-GeV) Au + Au (AGS) or Pb + Pb (SPS) colhsions in comparison 
to the experimental data from Refs. [53,23] for midrapidity (left column) and rapidity integrated 
yields (right column). The solid lines with open triangles show the results from the UrQMD 
calculations while the solid lines with open squares stem from the HSD approach with potentials. 
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FIG. 9. The excitation function of s and s production from central Au+Au (AGS) or Pb+Pb 
(SPS) collisions at impact parameter 6 = for HSD (l.h.s.) and UrQMD (r.h.s.) as appearing in 
mesons or baryons (antibaryons). The solid (dashed) lines with full (open) squares and circles 
indicate the number of s-quarks (s-quarks) in baryons (antibaryons) and mesons, respectively. 
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FIG. 10. The excitation function of the channel decomposition for and K~ from HSD 
(upper part) and UrQMD (lower part) from central Au + Au (AGS) or Ph + Ph (SPS) collisions 
at impact parameter 6 = 0. 



41 



p+p 



T ' 1 ' 1 ' 1 ' 1 ' 1 ' r 




FIG. 11. The rapidity distribution of protons from pp collisions at 4-160 GeV as calculated 
within HSD (solid lines) and UrQMD (dashed lines). The experimental data at 160 GeV are 
taken from Ref. [50]. 
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FIG. 12. The rapidity distribution of n^jK'^, K~ and A + S'^'s from pp collisions at 4-160 
GeV as calculated within HSD (solid lines) and UrQMD (dashed lines). The rapidity 
distributions from HSD and UrQMD at 4 and 8 GeV are scaled by factors of 6 and 3 while the 
K~ rapidity distributions are scaled by factors of 200 and 15, respectively. The experimental 
data for i^"*"'s and A + TP^s at 160 GeV are taken from Refs. [65] and [66]. 
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FIG. 13. The inclusive tt^ ,tt^ , , and A + S*^ production cross sections from pp colli- 
sions versus kinetic energy Eiab- The solid lines with open triangles show the UrQMD results, 
the solid line with full squares indicate the HSD results with 7^ defined by Eq. (2), whereas 
the dot-dashed lines correspond to the 7^ = 0.3. The exp. data are taken from Refs. [67] (full 
triangles), [68] (full and open circles), [69] (open squares) and [70] (stars). 
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FIG. 14. The rapidity distribution of ir^ , , K" and A + S'^'s from ir^p collisions at 2-8 
GeV/c as calculated within HSD (solid lines) and UrQMD (dashed lines). Note, that the 
and rapidity distributions from HSD and UrQMD at 2 GeV/c are scaled by factors of 5 and 
30, respectively. 
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FIG. 15. Total multiplicitis of 7r~^,i^+ and (i.e. 47r yields) from central Au + Au (at 
AGS) or P6 + Ph (at SPS) collisions in comparison to the total multiplicities from collisions 
(scaled by a factor 350/2) versus kinetic energy Eiah- The solid lines with full triangles and 
squares show the UrQMD (l.h.s.) and HSD results (r.h.s.) for AA collisions, respectively. The 
dotted lines with open triangles and squares correspond to the -pp multiplicities calculated within 
UrQMD (l.h.s.) and HSD (r.h.s.). 
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